Metagenomic Characterization and Volatile Compounds Determination in Rumen from Saanen Goat Kids Fed Olive Leaves.

The accumulation and disposal of by-products deriving from the agro-food industry represents a problem both from an economic and environmental point of view. The use of these matrices in zootechnical nutrition could represent a feasible solution. The aim of the study was to examine the effect of a diet containing olive leaves (OL), a by-product of the olive industry, on the ruminal microbial community of Saanen goat kids and on volatile organic compounds (VOCs) produced during the digestion. Twenty goat kids were randomly divided into two groups of ten goat kids each. The control group (CTR) was fed with a standard diet, while the experimental group (OL+) received a custom-formulated diet containing 10 % OL on a dry matter (DM) basis. After 30 days of trial, genomic DNA was extracted from the rumen liquor and prepared for 16S rRNA-gene sequencing to characterize the rumen microbiota; furthermore, rumen VOCs were also characterized by solid-phase microextraction coupled with gas chromatography-mass spectrometry. The Shannon's alpha index was not significantly different between the two groups, on the contrary, Bray-Curtis (p < 0.01) and Jaccard (p < 0.01) distances evidenced that feed affected microbial community. Eleven genera were influenced by OL supplementation, with a significant increase (p < 0.05) in Paludibacter, Fibrobacter, Sphaerochaeta Christensenella, Rikenella, Oligosphaera, Candidatus Endomicrobium, Anaerovorax, and Atopobium was observed, while the percentages of Bacteroides and Selenomonas were reduced (p < 0.05). Differences were also observed between the two groups at the family level (p < 0.004). Fibrobacteriaceae, Christensenellaceae, Coriobacteriaceae, Oligosphaeraceae, Candidatus Endomicrobium, and Planctomycetaceae were significantly higher (p < 0.05) in goat kids fed OL diet compared to CTR, while the levels of other identified families, Succinivibrionaceae and Bifidobacteriaceae, were opposite (p < 0.05). Finally, results showed that the main phyla in both groups were Bacteroidetes and Firmicutes; however, no significant differences in the relative abundance of any phyla were observed between the two groups. In addition to what has been reported, the analysis of VOCs at the rumen level showed the ability of the OL integration to induce an increase in hexanoic acid and a parallel decrease in decanal. Furthermore, only in OL+ samples there was the accumulation of α-terpineol to which a wide range of interesting biological properties is attributed. The presence of VOCs associated with health status suggests a favorable role of OL in preserving and improving animal welfare.

Metagenomic analysis reveals the shared and distinct features of the soil resistome across tundra, temperate prairie, and tropical ecosystems.

BACKGROUND: Soil is an important reservoir of antibiotic resistance genes (ARGs), but their potential risk in different ecosystems as well as response to anthropogenic land use change is unknown. We used a metagenomic approach and datasets with well-characterized metadata to investigate ARG types and amounts in soil DNA of three native ecosystems: Alaskan tundra, US Midwestern prairie, and Amazon rainforest, as well as the effect of conversion of the latter two to agriculture and pasture, respectively. RESULTS: High diversity (242 ARG subtypes) and abundance (0.184-0.242 ARG copies per 16S rRNA gene copy) were observed irrespective of ecosystem, with multidrug resistance and efflux pump the dominant class and mechanism. Ten regulatory genes were identified and they accounted for 13-35% of resistome abundances in soils, among them arlR, cpxR, ompR, vanR, and vanS were dominant and observed in all studied soils. We identified 55 non-regulatory ARGs shared by all 26 soil metagenomes of the three ecosystems, which accounted for more than 81% of non-regulatory resistome abundance. Proteobacteria, Firmicutes, and Actinobacteria were primary ARG hosts, 7 of 10 most abundant ARGs were found in all of them. No significant differences in both ARG diversity and abundance were observed between native prairie soil and adjacent long-term cultivated agriculture soil. We chose 12 clinically important ARGs to evaluate at the sequence level and found them to be distinct from those in human pathogens, and when assembled they were even more dissimilar. Significant correlation was found between bacterial community structure and resistome profile, suggesting that variance in resistome profile was mainly driven by the bacterial community composition. CONCLUSIONS: Our results identify candidate background ARGs (shared in all 26 soils), classify ARG hosts, quantify resistance classes, and provide quantitative and sequence information suggestive of very low risk but also revealing resistance gene variants that might emerge in the future. Video abstract.

Targeted assemblies of cas1 suggest CRISPR-Cas's response to soil warming.

There is an increasing interest in the clustered regularly interspaced short palindromic repeats CRISPR-associated protein (CRISPR-Cas) system to reveal potential virus-host dynamics. The universal and most conserved Cas protein, cas1 is an ideal marker to elucidate CRISPR-Cas ecology. We constructed eight Hidden Markov Models (HMMs) and assembled cas1 directly from metagenomes by a targeted-gene assembler, Xander, to improve detection capacity and resolve the diverse CRISPR-Cas systems. The eight HMMs were first validated by recovering all 17 cas1 subtypes from the simulated metagenome generated from 91 prokaryotic genomes across 11 phyla. We challenged the targeted method with 48 metagenomes from a tallgrass prairie in Central Oklahoma recovering 3394 cas1. Among those, 88 were near full length, 5 times more than in de-novo assemblies from the Oklahoma metagenomes. To validate the host assignment by cas1, the targeted-assembled cas1 was mapped to the de-novo assembled contigs. All the phylum assignments of those mapped contigs were assigned independent of CRISPR-Cas genes on the same contigs and consistent with the host taxonomies predicted by the mapped cas1. We then investigated whether 8 years of soil warming altered cas1 prevalence within the communities. A shift in microbial abundances was observed during the year with the biggest temperature differential (mean 4.16 °C above ambient). cas1 prevalence increased and even in the phyla with decreased microbial abundances over the next 3 years, suggesting increasing virus-host interactions in response to soil warming. This targeted method provides an alternative means to effectively mine cas1 from metagenomes and uncover the host communities.

Bioavailability of clay-adsorbed dioxin to Sphingomonas wittichii RW1 and its associated genome-wide shifts in gene expression.

Polychlorinated dibenzo-p-dioxins and dibenzofurans are a group of chemically-related pollutants categorically known as dioxins. Some of their chlorinated congeners are among the most hazardous pollutants that persist in the environment. This persistence is due in part to the limited number of bacteria capable of metabolizing these compounds, but also to their limited bioavailability in soil. We used Sphingomonas wittichii strain RW1 (RW1), one of the few strains able to grow on dioxin, to characterize its ability to respond to and degrade clay-bound dioxin. We found that RW1 grew on and completely degraded dibenzo-p-dioxin (DD) intercalated into the smectite clay saponite (SAP). To characterize the effects of DD sorption on RW1 gene expression, we compared transcriptomes of RW1 grown with either free crystalline DD or DD intercalated clay, i.e. sandwiched between the clay interlayers (DDSAP). Free crystalline DD appeared to cause greater expression of toxicity and stress related functions. Genes coding for heat shock proteins, chaperones, as well as genes involved in DNA repair, and efflux were up-regulated during growth on crystalline dioxin compared to growth on intercalated dioxin. In contrast, growth on intercalated dioxin up-regulated genes that might be important in recognition and uptake mechanisms, as well as surface interaction/attachment/biofilm formation such as extracellular solute-binding protein and LuxR. These differences in gene expression may reflect the underlying adaptive mechanisms by which RW1 cells sense and deploy pathways to access dioxin intercalated into clay. These data show that intercalated DD remains bioavailable to the degrading bacterium with implications for bioremediation alternatives.

Alterations in the Endophyte-Enriched Root-Associated Microbiome of Rice Receiving Growth-Promoting Treatments of Urea Fertilizer and Rhizobium Biofertilizer.

We examined the bacterial endophyte-enriched root-associated microbiome within rice (Oryza sativa) 55 days after growth in soil with and without urea fertilizer and/or biofertilization with a growth-promotive bacterial strain (Rhizobium leguminosarum bv. trifolii E11). After treatment to deplete rhizosphere/rhizoplane communities, washed roots were macerated and their endophyte-enriched communities were analyzed by 16S ribosomal DNA 454 amplicon pyrosequencing. This analysis clustered 99,990 valid sequence reads into 1105 operational taxonomic units (OTUs) with 97% sequence identity, 133 of which represented a consolidated core assemblage representing 12.04% of the fully detected OTU richness. Taxonomic affiliations indicated Proteobacteria as the most abundant phylum (especially α- and γ-Proteobacteria classes), followed by Firmicutes, Bacteroidetes, Verrucomicrobia, Actinobacteria, and several other phyla. Dominant genera included Rheinheimera, unclassified Rhodospirillaceae, Pseudomonas, Asticcacaulis, Sphingomonas, and Rhizobium. Several OTUs had close taxonomic affiliation to genera of diazotrophic rhizobacteria, including Rhizobium, unclassified Rhizobiales, Azospirillum, Azoarcus, unclassified Rhizobiaceae, Bradyrhizobium, Azonexus, Mesorhizobium, Devosia, Azovibrio, Azospira, Azomonas, and Azotobacter. The endophyte-enriched microbiome was restructured within roots receiving growth-promoting treatments. Compared to the untreated control, endophyte-enriched communities receiving urea and/or biofertilizer treatments were significantly reduced in OTU richness and relative read abundances. Several unique OTUs were enriched in each of the treatment communities. These alterations in structure of root-associated communities suggest dynamic interactions in the host plant microbiome, some of which may influence the well-documented positive synergistic impact of rhizobial biofertilizer inoculation plus low doses of urea-N fertilizer on growth promotion of rice, considered as one of the world's most important food crops.

Long-Term Effect of Different Fertilization and Cropping Systems on the Soil Antibiotic Resistome.

Different fertilization and cropping systems may influence short- and long-term residues of antibiotic resistance genes (ARGs) and mobile genetic elements (MGEs) in soil. Soils from dryland (peanut) and paddy (rice) fields, which originated from the same nonagricultural land (forested), were treated with either chemical fertilizer, composted manure, or no fertilizer for 26 years before sampling, which occurred one year after the last applications. ARGs and MGEs were investigated using highly parallel qPCR and high-throughput sequencing. Six of the 11 antibiotics measured by LC-MS/MS were detected in the manure applied soil, but not in the nonmanured soils, indicating their source was from previous manure applications. Compared to the unfertilized control, manure application did not show a large accumulation of ARGs in either cropping system but there were some minor effects of soil management on indigenous ARGs. Paddy soil showed higher accumulation of these ARGs, which corresponded to higher microbial biomass than the dryland soil. Chemical fertilizer increased relative abundance of these ARGs in dryland soil but decreased their relative abundance in paddy soil. These results show how long-term common soil management practices affect the abundance and type of ARGs and MGEs in two very different soil environments, one aerobic and the other primarily anaerobic.

Primer set 2.0 for highly parallel qPCR array targeting antibiotic resistance genes and mobile genetic elements.

The high-throughput antibiotic resistance gene (ARG) qPCR array, initially published in 2012, is increasingly used to quantify resistance and mobile determinants in environmental matrices. Continued utility of the array; however, necessitates improvements such as removing or redesigning questionable primer sets, updating targeted genes and coverage of available sequences. Towards this goal, a new primer design tool (EcoFunPrimer) was used to aid in identification of conserved regions of diverse genes. The total number of assays used for diverse genes was reduced from 91 old primer sets to 52 new primer sets, with only a 10% loss in sequence coverage. While the old and new array both contain 384 primer sets, a reduction in old primer sets permitted 147 additional ARGs and mobile genetic elements to be targeted. Results of validating the updated array with a mock community of strains resulted in over 98% of tested instances incurring true positive/negative calls. Common queries related to sensitivity, quantification and conventional data analysis (e.g. Ct cutoff value, and estimated genomic copies without standard curves) were also explored. A combined list of new and previously used primer sets is provided with a recommended set based on redesign of primer sets and results of validation.

ARGs-OAP v2.0 with an expanded SARG database and Hidden Markov Models for enhancement characterization and quantification of antibiotic resistance genes in environmental metagenomes.

Motivation: Much global attention has been paid to antibiotic resistance in monitoring its emergence, accumulation and dissemination. For rapid characterization and quantification of antibiotic resistance genes (ARGs) in metagenomic datasets, an online analysis pipeline, ARGs-OAP has been developed consisting of a database termed Structured Antibiotic Resistance Genes (the SARG) with a hierarchical structure (ARGs type-subtype-reference sequence). Results: The new release of the database, termed SARG version 2.0, contains sequences not only from CARD and ARDB databases, but also carefully selected and curated sequences from the latest protein collection of the NCBI-NR database, to keep up to date with the increasing number of ARG deposited sequences. SARG v2.0 has tripled the sequences of the first version and demonstrated improved coverage of ARGs detection in metagenomes from various environmental samples. In addition to annotation of high-throughput raw reads using a similarity search strategy, ARGs-OAP v2.0 now provides model-based identification of assembled sequences using SARGfam, a high-quality profile Hidden Markov Model (HMM), containing profiles of ARG subtypes. Additionally, ARGs-OAP v2.0 improves cell number quantification by using the average coverage of essential single copy marker genes, as an option in addition to the previous method based on the 16S rRNA gene. Availability and implementation: ARGs-OAP can be accessed through http://smile.hku.hk/SARGs. The database could be downloaded from the same site. Source codes for this study can be downloaded from https://github.com/xiaole99/ARGs-OAP-v2.0. Supplementary information: Supplementary data are available at Bioinformatics online.

Modulatory Influence of Segmented Filamentous Bacteria on Transcriptomic Response of Gnotobiotic Mice Exposed to TCDD.

Environmental toxicants such as 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), an aryl hydrocarbon receptor (AhR), are known to induce host toxicity and structural shifts in the gut microbiota. Key bacterial populations with similar or opposing functional responses to AhR ligand exposure may potentially help regulate expression of genes associated with immune dysfunction. To examine this question and the mechanisms for AhR ligand-induced bacterial shifts, C57BL/6 gnotobiotic mice were colonized with and without segmented filamentous bacteria (SFB) - an immune activator. Mice were also colonized with polysaccharide A producing Bacteroides fragilis - an immune suppressor to serve as a commensal background. Following colonization, mice were administered TCDD (30 µg/kg) every 4 days for 28 days by oral gavage. Quantified with the nCounter® mouse immunology panel, opposing responses in ileal gene expression (e.g., genes associated with T-cell differentiation via the class II major histocompatibility complex) as a result of TCDD dosing and SFB colonization were observed. Genes that responded to TCDD in the presence of SFB did not show a significant response in the absence of SFB, and vice versa. Regulatory T-cells examined in the mesenteric lymph-nodes, spleen, and blood were also less impacted by TCDD in mice colonized with SFB. TCDD-induced shifts in abundance of SFB and B. fragilis compared with previous studies in mice with a traditional gut microbiome. With regard to the mouse model colonized with individual populations, results indicate that TCDD-induced host response was significantly modulated by the presence of SFB in the gut microbiome, providing insight into therapeutic potential between AhR ligands and key commensals.

Isothermal assay targeting class 1 integrase gene for environmental surveillance of antibiotic resistance markers.

Antimicrobial resistance genes (ARGs) present in the environment pose a risk to human health due to potential for transfer to human pathogens. Surveillance is an integral part of mitigating environmental dissemination. Quantification of the mobile genetic element class 1 integron-integrase gene (intI1) has been proposed as a surrogate to measuring multiple ARGs. Measurement of such indicator genes can be further simplified by adopting emerging nucleic acids methods such as loop mediated isothermal amplification (LAMP). In this study, LAMP assays were designed and tested for estimating relative abundance of the intI1 gene, which included design of a universal bacteria 16S rRNA gene assay. Following validation of sensitivity and specificity with known bacterial strains, the assays were tested using DNA extracted from river and lake samples. Results showed a significant Pearson correlation (R2 = 0.8) between the intI1 gene LAMP assay and ARG relative abundance (measured via qPCR). To demonstrate the ruggedness of the LAMP assays, experiments were also run in the hands of relatively "untrained" personnel by volunteer undergraduate students at a local community college using a hand-held real-time DNA analysis device - Gene-Z. Overall, results support use of the intI1 gene as an indicator of ARGs and the LAMP assays exhibit the opportunity for volunteers to monitor environmental samples for anthropogenic pollution outside of a specialized laboratory.

Virulence factor activity relationships (VFARs): a bioinformatics perspective.

Virulence factor activity relationships (VFARs) - a concept loosely based on quantitative structure-activity relationships (QSARs) for chemicals was proposed as a predictive tool for ranking risks due to microorganisms relevant to water safety. A rapid increase in sequencing capabilities and bioinformatics tools has significantly increased the potential for VFAR-based analyses. This review summarizes more than 20 bioinformatics databases and tools, developed over the last decade, along with their virulence and antimicrobial resistance prediction capabilities. With the number of bacterial whole genome sequences exceeding 241 000 and metagenomic analysis projects exceeding 13 000 and the ability to add additional genome sequences for few hundred dollars, it is evident that further development of VFARs is not limited by the availability of information at least at the genomic level. However, additional information related to co-occurrence, treatment response, modulation of virulence due to environmental and other factors, and economic impact must be gathered and incorporated in a manner that also addresses the associated uncertainties. Of the bioinformatics tools, a majority are either designed exclusively for virulence/resistance determination or equipped with a dedicated module. The remaining have the potential to be employed for evaluating virulence. This review focusing broadly on omics technologies and tools supports the notion that these tools are now sufficiently developed to allow the application of VFAR approaches combined with additional engineering and economic analyses to rank and prioritize organisms important to a given niche. Knowledge gaps do exist but can be filled with focused experimental and theoretical analyses that were unimaginable a decade ago. Further developments should consider the integration of the measurement of activity, risk, and uncertainty to improve the current capabilities.

RNAcentral: a comprehensive database of non-coding RNA sequences.

RNAcentral is a database of non-coding RNA (ncRNA) sequences that aggregates data from specialised ncRNA resources and provides a single entry point for accessing ncRNA sequences of all ncRNA types from all organisms. Since its launch in 2014, RNAcentral has integrated twelve new resources, taking the total number of collaborating database to 22, and began importing new types of data, such as modified nucleotides from MODOMICS and PDB. We created new species-specific identifiers that refer to unique RNA sequences within a context of single species. The website has been subject to continuous improvements focusing on text and sequence similarity searches as well as genome browsing functionality. All RNAcentral data is provided for free and is available for browsing, bulk downloads, and programmatic access at http://rnacentral.org/.

Influence of Soil Characteristics and Proximity to Antarctic Research Stations on Abundance of Antibiotic Resistance Genes in Soils.

Soil is an important environmental reservoir of antibiotic resistance genes (ARGs), which are increasingly recognized as environmental contaminants. Methods to assess the risks associated with the acquisition or transfer of resistance mechanisms are still underdeveloped. Quantification of background levels of antibiotic resistance genes and what alters those is a first step in understanding our environmental resistome. Toward this goal, 62 samples were collected over 3 years from soils near the 30-year old Gondwana Research Station and for 4 years before and during development of the new Jang Bogo Research Station, both at Terra Nova Bay in Antarctica. These sites reflect limited and more extensive human impact, respectively. A qPCR array with 384 primer sets targeting antibiotic resistance genes and mobile genetic elements (MGEs) was used to detect and quantify these genes. A total of 73 ARGs and MGEs encompassing eight major antibiotic resistance gene categories were detected, but most at very low levels. Antarctic soil appeared to be a common reservoir for seven ARGs since they were present in most samples (42%-88%). If the seven widespread genes were removed, there was a correlation between the relative abundance of MGEs and ARGs, more typical of contaminated sites. There was a relationship between ARG content and distance from both research stations, with a significant effect at the Jang Bogo Station especially when excluding the seven widespread genes; however, the relative abundance of ARGs did not increase over the 4 year period. Silt, clay, total organic carbon, and SiO2 were the top edaphic factors that correlated with ARG abundance. Overall, this study identifies that human activity and certain soil characteristics correlate with antibiotic resistance genes in these oligotrophic Antarctic soils and provides a baseline of ARGs and MGEs for future comparisons.

Sequence-based classification and identification of Fungi.

Fungal taxonomy and ecology have been revolutionized by the application of molecular methods and both have increasing connections to genomics and functional biology. However, data streams from traditional specimen- and culture-based systematics are not yet fully integrated with those from metagenomic and metatranscriptomic studies, which limits understanding of the taxonomic diversity and metabolic properties of fungal communities. This article reviews current resources, needs, and opportunities for sequence-based classification and identification (SBCI) in fungi as well as related efforts in prokaryotes. To realize the full potential of fungal SBCI it will be necessary to make advances in multiple areas. Improvements in sequencing methods, including long-read and single-cell technologies, will empower fungal molecular ecologists to look beyond ITS and current shotgun metagenomics approaches. Data quality and accessibility will be enhanced by attention to data and metadata standards and rigorous enforcement of policies for deposition of data and workflows. Taxonomic communities will need to develop best practices for molecular characterization in their focal clades, while also contributing to globally useful datasets including ITS. Changes to nomenclatural rules are needed to enable validPUBLICation of sequence-based taxon descriptions. Finally, cultural shifts are necessary to promote adoption of SBCI and to accord professional credit to individuals who contribute to community resources.

Transcript Profile of Flowering Regulatory Genes in VcFT-Overexpressing Blueberry Plants.

In order to identify genetic components in flowering pathways of highbush blueberry (Vaccinium corymbosum L.), a transcriptome reference composed of 254,396 transcripts and 179,853 gene contigs was developed by assembly of 72.7 million reads using Trinity. Using this transcriptome reference and a query of flowering pathway genes of herbaceous plants, we identified potential flowering pathway genes/transcripts of blueberry. Transcriptome analysis of flowering pathway genes was then conducted on leaf tissue samples of transgenic blueberry cv. Aurora ('VcFT-Aurora'), which overexpresses a blueberry FLOWERING LOCUS T-like gene (VcFT). Sixty-one blueberry transcripts of 40 genes showed high similarities to 33 known flowering-related genes of herbaceous plants, of which 17 down-regulated and 16 up-regulated genes were identified in 'VcFT-Aurora'. All down-regulated genes encoded transcription factors/enzymes upstream in the signaling pathway containing VcFT. A blueberry CONSTANS-LIKE 5-like (VcCOL5) gene was down-regulated and associated with five other differentially expressed (DE) genes in the photoperiod-mediated flowering pathway. Three down-regulated genes, i.e., a MADS-AFFECTING FLOWERING 2-like gene (VcMAF2), a MADS-AFFECTING FLOWERING 5-like gene (VcMAF5), and a VERNALIZATION1-like gene (VcVRN1), may function as integrators in place of FLOWERING LOCUS C (FLC) in the vernalization pathway. Because no CONSTAN1-like or FLOWERING LOCUS C-like genes were found in blueberry, VcCOL5 and VcMAF2/VcMAF5 or VRN1 might be the major integrator(s) in the photoperiod- and vernalization-mediated flowering pathway, respectively. The major down-stream genes of VcFT, i.e., SUPPRESSOR of Overexpression of Constans 1-like (VcSOC1), LEAFY-like (VcLFY), APETALA1-like (VcAP1), CAULIFLOWER 1-like (VcCAL1), and FRUITFULL-like (VcFUL) genes were present and showed high similarity to their orthologues in herbaceous plants. Moreover, overexpression of VcFT promoted expression of all of these VcFT downstream genes. These results suggest that VcFT's down-stream genes appear conserved in blueberry.

Sphingomonas wittichii Strain RW1 Genome-Wide Gene Expression Shifts in Response to Dioxins and Clay.

Sphingomonas wittichii strain RW1 (RW1) is one of the few strains that can grow on dibenzo-p-dioxin (DD). We conducted a transcriptomic study of RW1 using RNA-Seq to outline transcriptional responses to DD, dibenzofuran (DF), and the smectite clay mineral saponite with succinate as carbon source. The ability to grow on DD is rare compared to growth on the chemically similar DF even though the same initial dioxygenase may be involved in oxidation of both substrates. Therefore, we hypothesized the reason for this lies beyond catabolic pathways and may concern genes involved in processes for cell-substrate interactions such as substrate recognition, transport, and detoxification. Compared to succinate (SUC) as control carbon source, DF caused over 240 protein-coding genes to be differentially expressed, whereas more than 300 were differentially expressed with DD. Stress response genes were up-regulated in response to both DD and DF. This effect was stronger with DD than DF, suggesting a higher toxicity of DD compared to DF. Both DD and DF caused changes in expression of genes involved in active cross-membrane transport such as TonB-dependent receptor proteins, but the patterns of change differed between the two substrates. Multiple transcription factor genes also displayed expression patterns distinct to DD and DF growth. DD and DF induced the catechol ortho- and the salicylate/gentisate pathways, respectively. Both DD and DF induced the shared down-stream aliphatic intermediate compound pathway. Clay caused category-wide down-regulation of genes for cell motility and chemotaxis, particularly those involved in the synthesis, assembly and functioning of flagella. This is an environmentally important finding because clay is a major component of soil microbes' microenvironment influencing local chemistry and may serve as a geosorbent for toxic pollutants. Similar to clay, DD and DF also affected expression of genes involved in motility and chemotaxis.

ARGs-OAP: online analysis pipeline for antibiotic resistance genes detection from metagenomic data using an integrated structured ARG-database.

MOTIVATION: Environmental dissemination of antibiotic resistance genes (ARGs) has become an increasing concern for public health. Metagenomics approaches can effectively detect broad profiles of ARGs in environmental samples; however, the detection and subsequent classification of ARG-like sequences are time consuming and have been severe obstacles in employing metagenomic methods. We sought to accelerate quantification of ARGs in metagenomic data from environmental samples. RESULTS: A Structured ARG reference database (SARG) was constructed by integrating ARDB and CARD, the two most commonly used databases. SARG was curated to remove redundant sequences and optimized to facilitate query sequence identification by similarity. A database with a hierarchical structure (type-subtype-reference sequence) was then constructed to facilitate classification (assigning ARG-like sequence to type, subtype and reference sequence) of sequences identified through similarity search. Utilizing SARG and a previously proposed hybrid functional gene annotation pipeline, we developed an online pipeline called ARGs-OAP for fast annotation and classification of ARG-like sequences from metagenomic data. We also evaluated and proposed a set of criteria important for efficiently conducting metagenomic analysis of ARGs using ARGs-OAP. AVAILABILITY AND IMPLEMENTATION: Perl script for ARGs-OAP can be downloaded from https://github.com/biofuture/Ublastx_stageone ARGs-OAP can be accessed through http://smile.hku.hk/SARGs CONTACT: zhangt@hku.hk or tiedjej@msu.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

RNAcentral: an international database of ncRNA sequences.

The field of non-coding RNA biology has been hampered by the lack of availability of a comprehensive, up-to-date collection of accessioned RNA sequences. Here we present the first release of RNAcentral, a database that collates and integrates information from an international consortium of established RNA sequence databases. The initial release contains over 8.1 million sequences, including representatives of all major functional classes. A web portal (http://rnacentral.org) provides free access to data, search functionality, cross-references, source code and an integrated genome browser for selected species.

Gene expression analysis of E. coli strains provides insights into the role of gene regulation in diversification.

Escherichia coli spans a genetic continuum from enteric strains to several phylogenetically distinct, atypical lineages that are rare in humans, but more common in extra-intestinal environments. To investigate the link between gene regulation, phylogeny and diversification in this species, we analyzed global gene expression profiles of four strains representing distinct evolutionary lineages, including a well-studied laboratory strain, a typical commensal (enteric) strain and two environmental strains. RNA-Seq was employed to compare the whole transcriptomes of strains grown under batch, chemostat and starvation conditions. Highly differentially expressed genes showed a significantly lower nucleotide sequence identity compared with other genes, indicating that gene regulation and coding sequence conservation are directly connected. Overall, distances between the strains based on gene expression profiles were largely dependent on the culture condition and did not reflect phylogenetic relatedness. Expression differences of commonly shared genes (all four strains) and E. coli core genes were consistently smaller between strains characterized by more similar primary habitats. For instance, environmental strains exhibited increased expression of stress defense genes under carbon-limited growth and entered a more pronounced survival-like phenotype during starvation compared with other strains, which stayed more alert for substrate scavenging and catabolism during no-growth conditions. Since those environmental strains show similar genetic distance to each other and to the other two strains, these findings cannot be simply attributed to genetic relatedness but suggest physiological adaptations. Our study provides new insights into ecologically relevant gene-expression and underscores the role of (differential) gene regulation for the diversification of the model bacterial species.

Ribosomal Database Project: data and tools for high throughput rRNA analysis.

Ribosomal Database Project (RDP; http://rdp.cme.msu.edu/) provides the research community with aligned and annotated rRNA gene sequence data, along with tools to allow researchers to analyze their own rRNA gene sequences in the RDP framework. RDP data and tools are utilized in fields as diverse as human health, microbial ecology, environmental microbiology, nucleic acid chemistry, taxonomy and phylogenetics. In addition to aligned and annotated collections of bacterial and archaeal small subunit rRNA genes, RDP now includes a collection of fungal large subunit rRNA genes. RDP tools, including Classifier and Aligner, have been updated to work with this new fungal collection. The use of high-throughput sequencing to characterize environmental microbial populations has exploded in the past several years, and as sequence technologies have improved, the sizes of environmental datasets have increased. With release 11, RDP is providing an expanded set of tools to facilitate analysis of high-throughput data, including both single-stranded and paired-end reads. In addition, most tools are now available as open source packages for download and local use by researchers with high-volume needs or who would like to develop custom analysis pipelines.